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Abstract 

A reexamination of hadronic comover scattering indicates that this mechanism 
cannot explain the observed ip suppression in Pb+Pb interactions. The possibihty 
of quark-gluon plasma formation is therefore considered. 



1 Introduction 

Ever since it was realized that the ip suppression by quark-gluon plasma 
formation predicted by Matsui and Satz [Q] could be mimicked by hadronic 
means, the source of the observed if) suppression has remained controver- 
sial p, ^, ^, ^, ^. In this paper we address the compatibility of the data 
with hadronic models in light of a reanalysis of our previous work We 
also discuss a possible quark-gluon plasma origin of the apparent suppression 
thresholds in the latest Pb-fPb data [jlO . 
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Department of Energy under Contract No. DE-AC03-76SF0098. 
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2 A Dependence of Hadronic Models 



In this section, we show that the most simphstic assumptions of ip suppression 
in hadronic interactions lead to the same A dependence for absorption and 
secondary, or comover, scatterings. In hadron-nucleus, hA, collisions [jlT| 



where b is the impact parameter, T| (&) is the effective nuclear profile 



Tfib) 



dz pA{b, z) exp 



dz'pA^b, z')a^N{z' - z] 



(2) 



and a^N is the cross section for ip (or \ccg)) absorption by nucleons. The 
nuclear density distributions, pA{bi z), are taken from nuclear charge density 
measurements |T^. The comover survival probability is |Tl[ 



exp 



/ dr {a^cov)n{T, b) 

J Tn 



exp 



(3) 



where a^co is the -i/^-comover absorption cross section, v ~ 0.6 is the relative 
■j/'-comover velocity, and n(r, b) is the density of comovers at time r and 
impact parameter b. Integrating over r and relating the initial density to the 
final hadron rapidity density JTSf, nQT^ = {it Rl)~^{dN/dy), we have [p] 



{(^iicov) / dTn{T,b) 



dN. 



2 (ThNTA{b) — \y^^ 



ln(^ 

To 



(4) 



where dN/dy\y^^, the central rapidity density in an hp collision, is scaled up to 
hA interactions by ahNTA^b), the number of participants. Since nRl ?s ahN, 
these two factors cancel. The effective proper lifetime tj over which the 
comovers, formed at time tq, interact with the is tj ^ Rh/v. 

Expanding the exponentials in eqs. @ and @ to terms linear in cross 
section, integrating eq. (|1|) and re-exponentiating, we obtain |TT 



Aa 



exp 



IQlTTn 



, , , fTi\ dN . 
2(w)ln(^-j— I,, 



exp 



{-A^'^{^ + (3)] . 



(5) 
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For large targets with A > 50, A^^^ ^ In A so that 



(7) 



The comover contribution to a could be significant |]TT| unless a^co is strongly 
reduced at low energies [1^ (although see also JlSl), tj is very short or dN/dy 
is small. In most recent analyses p, ^, ||, |I6|, an effective value of a^N 
including both effects is determined from a. 

Unfortunately the identical nature of the two contributions to eq. (|^) sug- 
gests that they may be inextricably intertwined. Experimental constraints 
on cr^Tv and a^co would clearly be valuable. Exclusive ip production in near 
threshold pA interactions fl^ and an inverse kinematics experiment [|18| 
could clarify ct^at. A measurement of dN/dy associated with ip production 
could place limits on the comover contribution. 



3 Et Dependence of Hadronic Models 

The transverse energy, E^, dependence of production in nucleus- nucleus 
collisions is generally compared to Drell-Yan, DY, production, uninfluenced 
by flnal-state interactions. The ip /DY ratio as a function of is 

Ba^ ^^^^ ^ Ba,,^^Jd'bJ dhTf jb) Tf{\b - s\) SABjb, s) p{ET;b) 



' app^^^^- J d%J d^sTAib)TBi\b- s\)piEr,b) 

The probability to produce Et at impact parameter b, p{Et] b), is a Gaussian 
with mean Exib) = eNN^Bib) proportional to the number of participants and 
standard deviation cr^(6) = uje]^ET{b) [llj. The parameters e^-, the energy 
per participant, and uo are chosen to agree with the NA38/NA50 neutral Et 
distributions. In S+U interactions with 1.7 < r] < 4.1 |]19|, eN = 0.74 GeV 
while in Pb-^Pb interactions with 1.1 < r] < 2.3 [H, eN = 0.4 GeV. We use 
uj = 3.2 but note that the precise value of uj only weakly influences the ip/DY 
ratios. In Fig. |l| the Drell-Yan Et distributions are shown normalized to the 
pp cross section at 200 GeV. (See also pO| .) 

The ijj, ijj' and Xc are assumed to interact with nucleons while in \ccg) color 
octet states with a lifetime of ~ 0.3 fm [^. Because the flnal charmonium 



state has formed by the time it interacts with comovers, the ~ 30% Xc 
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Figure 1: The Et distributions of Drell-Yan pairs with 2.9 < M < 4.5 GeV 
in (a) S+U and (b) Pb+Pb interactions. The pp normahzation is calculated 
at 200 GeV with the GRV LO 0] distributions and I cos^csl < 0.5. 



and ~ 12% ip' decay contributions to ijj production are considered 
separately. To include V'' and Xc feeddown to the ip, we scale the ip' and Xc 
comover cross sections relative to a^co by the squares of the radii ||25| so that 
o-^'co ~ 3.8a^co, c^xcco ~ 2.4cr^co and 

SAB{b, s) = 0.58 S^b, s) + 0.3 S^^{b, s) + 0.12 S^>{b, s) . (9) 

In Refs. 0, |l^ the density of comovers was assumed to be directly propor- 
tional to Et, tlco = uEt/ Et{0) with tq ~ 2 fm and tj ^ Ra/v. However, 
since Et is proportional to the number of participants, we replace hoTq by 
nAB{b,s), the participant density 0, 

SAB{b, s) = expi^-{acov)nAB{b, s) In | , (10) 

similar to Ref . . The comover formation time, a function of the path length 
L H, 13, ro(6, s) = 1 + La(&, s)/7a(&, s) + s)/7b(6, s), is ~ 2 fm in 
central collisions and ~ 1 fm in very peripheral collisions. The comovers 
interact with the ip only if t/(6) > ro(6, s) where 

(U\^\^a/v b<RB-RA /..N 

\ {RA + RB-b)/{2v) RB-RA<b<RB + RA ■ ^ ' 
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Figure ^ shows the ip/DY ratios for ct^n = 7.3 mb and 4.8 mb. Comovers 
are included when a^^ = 4.8 mb. Previously we used a^co ~ 2cr^Ar/3 = 



3.2 mb from quark-counting arguments; here we find a^co = 0.67 mb fits the 
S+U data, similar to the results of Ref. 0. When feeddown is not included, 
a^co = 1-1 nib produces equivalent agreement. The ip cross section in pp 
interactions is fit to the A dependence given each value of a^^: (Jpp^^ = 2.1 
nb for (J^N = 4.8 mb and (Jpp^^ = 2.3 nb for c^at = 7.3 mb. Assuming 
the maximum comover pA contribution when a^N = 4.8 mb is equivalent to 
an effective a^N, eqs. of ~ 6 mb, in agreement with the pA data [@]. 



The values of app^^ are in agreement with both the data and low energy 



parameterizations of the cross section 0, |2^ as well as with next-to-leading 
order calculations pi|. 
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Figure 2: The tp/DY calculations are compared with data [Q] from S-l-U 
(a) and Pb-|-Pb (b) interactions, both at 200 GeV. The solid curve shows 
a^N = 4.8 mb and a^co = 0.67 mb. The dashed and dot-dashed assume 
o'l/ico = for a^N = 4.8 mb and a^N = 7.3 mb, respectively. 



The results in Fig. ^ differ from those obtained in Ref. where simul- 
taneous agreement of the model with both the S-l-U and Pb-|-Pb data was 
found. The agreement with the NA38 S-l-U data H has improved over that 
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of Refs. 1^, |T6[. However, the result now disagrees with the Pb+Pb data, 
in accord with the conclusions of Refs. 0, The major difference lies in 
the normalization of the ip/DY ratio, determined from the pp production 
cross sections in the NA50 phase space. All the data [0, has been 

isospin adjusted to the pp cross section, continuum mass to 2.9 < M < 4.5 
GeV, projectile energy to 200 GeV, and phase space to < j/cm < 1 and 
I cos 6*03 1 < 0.5. The leading order Drell-Yan cross sections were adjusted 
for mass and isospin using the GRV LO distributions with the NA50 K 
factor 0. The isospin correction at = 19.4 GeV and 2.9 < M < 4.5 
GeV is pp/Pb+Pb ~ 1.3. Since ip production is dominated by gluon fusion, 
it is insensitive to isospin. In Refs. 0, |lT], |16| the S+U data was compared 
to the "continuum" mass range, 1.7 < M < 2.7 GeV, and extrapolated to 
I coshes I < 1- Contrary to what was stated in Ref. 0], the angular ad- 
justment from |cos6'cs| < 1 to | coshes I < 0-5 was left out of the Pb+Pb 
calculation. Thus a 23% increase of the normalization in Ref. H is needed 



because while the ip decays decays isotropically to leptons requiring a 
factor of two adjustment, the Drell-Yan cross section, oc 1 + cos^6'cs [0, is 
adjusted by 2.46. Therefore, although assuming a^N = 4.8 mb still leads to 
the earlier conclusion that comovers are necessary to explain the S+U data 



[ p!T| , |T6[ , the Pb+Pb results now suggest that the more absorption is needed. 



Other recent comover calculations |^ find simultaneous agreement be- 
tween the two systems within dynamical models of secondary production, 
consistent with scaling npb+pb by a factor of two. If secondary production 
in the central rapidity region increased oc where x > 1, as suggested by 
recent results from NA49 the discrepancy between our results and the 



Pb+Pb data could be reduced at high Et, e.g. if a; = 1.3, ■j/'/DY decreases 
13% at Et = 150 GeV. Agreement of the comover interpretation with the 
data was also found in Ref. [Q. While some of the differences with our con- 
clusions lie in the choice of model parameters and the specific comover Ej^ 
dependence, the results are also rather sensitive to the nuclear density pro- 
file. The S+U pT-dependent data favor a Woods-Saxon density profile over 
a sharp surface nucleus However, the precise Woods-Saxon may also 

affect the agreement of the model with the data. For example, Ref. uses 
a profile with lower central nuclear densities |33] than those determined from 
data |T^. In particular, ps(0,0) is ~ 35% lower while pu(0,0) and /9pb(0,0) 
10% lower. 

The ip'/ip ratios are shown in Fig. ^. If a^'co is a free parameter. 
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Figure 3: The calculated ip'/ip ratio is compared with data from S+U (a) 
and Pb+Pb (b) collisions |3^. The solid curves show a^N = ct^'n = 4.8 
mb, a^co = 0.67 mb and a^'co = 3.8cr^co- The dashed curves, neglecting 
feeddown, take a^'co = 8.5 mb. 



neglecting feeddown, then cr^'co ~ 7-8.5 mb reproduces the observed sup- 
pression for a^N = 7.3 mb with a^co ~ and a^N = 4.8 mb, a^co = 0.67 
mb respectively. Including feeddown, eq. @, a^'co a^co and the calculated 
ratios underpredict the data. Other combinations of cr^co, c"xcco and a^'co 
produce similar -^/DY results but do not significantly improve agreement 
with the if)' / ip ratios. 



4 Suppression by Plasma Screening? 

In light of these hadronic results, we speculate as to the nature of ip suppres- 
sion by plasma screening. The quarkonium potential is expected to be mod- 
ified at finite temperatures by the screening mass /i(T). When /i(T) > fiD, 
screening prevents resonance formation at T = T^. The critical values are 
1^1' = 0.342 GeV, /i^' = 0.357 GeV and /x^ = 0.699 GeV [0. Whether a 
simple plasma interpretation is possible depends strongly on /x(T). Plasma 
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suppression is possible when the energy density, e, is greater than en- 

We assume that /i(T) = Jl + nf/6g{T)T where g^{T) = 487rV[(33 - 



2nf)\nF^] with F = K{T/T^){TjAj^) |3|. In SU(3) gauge theory, = 
260 MeV and TJA^ = 1.03 ± 0.19 |8[. A fit to the heavy quark 



potential at high temperatures yields K ^ 33.8 |^. Lattice results with 
Uf =2 and 4 suggest T, = 170 MeV and TJA^ = 1.05 If the SU(3) 
value of K is applicable when Uf > 0, then for = 2 or 3, the Xc breaks up 
at T^" ^ 180-200 MeV and eg= = 2.1-2.6 GeV/fm^ while the ijj' subsequently 
breaks up at T^' ^ 190-210 MeV and e^' = 2.6-3.1 GeV/fm^. The ^ itself 
would not break up until T > 2Tc. If = 4, = Tc for both the Xc and 

with = e, ^ 2.1 GeV/fm^ while ^ 370 MeV. The Uf = results 
are similar to those with nj = A albeit with correspondingly larger e due to 
the higher Tf.. However, the high temperature limit is probably invalid for 
T/Tc < 3.5 |3^. A fit of to lattice data for T > Tc results in larger 
values of /^(T), suggesting To = Tc for the Xc ^' and ip with e£) = ec ~ 1.3-2.1 
GeV/fm^ for nj = 2-4 and = Cc ~ 3.1 GeV/fm^ for ri/ = 0. 

Thus the ip/DY ratio could contain one or two thresholds, as suggested 
by the new Pb-|-Pb data [|1^. We choose three illustrative cases: I) Uf = 3, 
sequential Xc and ip' break up; II) nj = 4, T^" = = T^, and III) Tq = 
for all charmonium states. Cases I and II assume K ^ 33.8 while case III 
takes K from the fit for T > Q. 

The relationship between e, not directly measurable, and Et is deter- 
mined from n{ET)i the average number of subcollisions per area 0, p^] . 



Using the NA50 prescription |^ to calculate e in Pb-|-Pb collisions places 
the maximum e in S-l-U collisions below the first threshold in the Pb-f-Pb 
data 0. 

If suppression is assumed to be total at e = en, then sharp thresholds 
appear, as shown for case I in Fig. §(a). Such sharp thresholds seem physi- 
cally unlikely due to the fluctuations in Et-, the finite size of the plasma, and 
the finite momentum of the charmonium state. The smooth increase of the 
suppression for e > e^) in a plasma with radius i?pb and {) is shown in 
the remaining curves. All calculations assume a^N = 4.8 mb and a^co = 0.67 
mb with eq. (P) modified to include plasma suppression, 

SABib, s) = 0.58 ^J^P(e) + 0.3 Sjf ^(e) + 0.12 S^P^'ie) . (12) 

A sudden change of slope in /DY appears when the plasma suppression 
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Figure 4: The Pb+Pb ip/DY § (a) and ip'/ip Q (b) ratios at 200 GeV 
compared to plasma predictions assuming a^N = ct^'n = 4.8 mb, a^co = 0.67 
mb and a^'co = 3.8cr^co- The sohd and dashed curves illustrate case I with 
total suppression above en and with R = i?pb respectively. The dot-dashed 
and dotted curves correspond to cases II and III with R = i?pb. 



begins, even without assuming a sharp threshold, not predicted by hadronic 
models. Note that case III overpredicts the suppression and that the thresh- 



olds introduced in case I are somewhat low relative to |T0[ . Increasing px with 
R = i?pb does not significantly change the suppression pattern. If i? = 1 fm, 
the plasma contribution is negligible for > 4 GeV in case I and II while 
case III is comparable to the dashed and dot-dashed curves when p-r ~ 3 
GeV. Since (pt) ~ 1 GeV at y/s = 19.4 GeV the magnitude of the sup- 
pression suggests that only the Xc and ip' are suppressed at current energies. 
Similar results are obtained if a^Tv = 7.3 mb and comovers are not included. 

The ip'/i/j ratio assuming feeddown is shown in Fig. ^(b). If a sharp 
threshold occurs, such as in case I, the ip'/ip ratio vanishes above and 
abruptly increases at , an unlikely scenario in a real collision. Case II 
produces agreement with the data but cannot explain why the feeddown 
scenario disagrees with the S+U data. An enhanced excitation of ip into ip' 
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near the chiral transition |^ could also decrease the ip' /ip ratio. 



5 Summary 

We have reviewed ip suppression by hadronic means and discussed the Et 
dependence of a plasma component in Pb+Pb interactions. We showed that 
the naive A dependencies of comover scattering and nuclear absorption are 
identical. We have extended our comover calculation and found that the 
S+U and Pb+Pb data are not simultaneously described in our approach 
although agreement in other scenarios with different models of secondary 
production cannot be ruled out ^. Thus ip suppression by plasma 

production should also be considered. 

The plasma models depend strongly on /i(T). Sharp thresholds do not 
occur unless introduced, as shown also in ||^, In particular, assuming a 
central plasma density proportional to UAsip, s) 0, only an increase in slope 
occurs at the onset of the plasma phase. The Et and e correlation needs to be 
clarified, especially since Et fluctuations cause the e bins to overlap [l^. A 



positive plasma indication also needs to be confirmed in other measurements. 
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